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The precise data for the total cross section a(e'*'e~ hadrons) from the charm threshold region, when combined 
with the evaluation of moments with three loop accuracy, lead to a direct determination of the short distance 
MS charm quark mass mc(mc) = 1.304(27). Applying the same approach to the bottom quark we obtain 
mb{mh) = 4.191(51) GeV. A complementary method for the determination of mb is based of the analysis of the 
T(IS') system which is confronted with a next-to-next-to-next-to-leading order calculation of the corresponding 
energy level. This leads to mb(mb) = 4.346(70) GeV. 


1. Introduction 

During the past years new and more precise 
data for CT(e“'"e“ —> hadrons) have become avail¬ 
able in the low energy region between 2 and 
10 GeV. At the same time increasingly precise cal¬ 
culations have been performed in the framework 
of perturbative QGD (pQGD), both for the cross 
section as a function of the center-of-mass energy 
t/s, including quark mass effects, and for its mo¬ 
ments which allow for a precise determination of 
the quark mass. A fresh look at the evaluation of 
the charm quark mass with the help of sum rules 
is thus an obvious task. We will concentrate on 
low moments as suggested by the ITEP group 
long ago [ Q. This is a natural route to deter¬ 
mine directly a short distance mass, say mc(TOc), 
in the MS scheme [ ||. 

The same method can also be applied to the 
bottom system leading directly to mb{mb). 

A complementary method for the determina¬ 
tion of mb is based of the analysis of the T(15') 
system. Recently a major improvement on the 
theoretical side has been achieved by the evalua¬ 
tion of the next-to-next-to-next-to-leading order 
(N^LO) corrections to the energy level. The dis¬ 
advantage of this method is due to the large non- 
perturbative effects. 
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2002 . 


2. Charm and bottom quark mass from 
low-order moment sum rules 


A very elegant method for the determination of 
the charm quark mass is based on the direct com¬ 
parison of theoretical and experimental moments 
of the charm quark contribution to the photon 
polarization function. In the limit of small mo¬ 
mentum the latter can be cast into the form [ 


nc(g') 



n>0 


( 1 ) 


with z = q'^Kiml) where rric = md/J-) is the MS 
charm quark mass at the scale /r. The perturba¬ 
tive series for the coefficients Cji up to n = 8 is 
known analytically [ ^ up to order a^. We define 
the moments through 




Md = IqI 



( 2 ) 


which leads to the following formula for the charm 
quark mass 


md^^) 


1 / C„(lnme) \ 

2 V J 


( 3 ) 


The analysis of the experimental moments with 
n = 1,...,4 leads to the results displayed in 

Fig.0[|l. 
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The moment with n = 1 is evidently least sensi¬ 
tive to non-perturbative contributions from con¬ 
densates, to the Coulombic higher order effects, 
the variation of /r and the parametric as depen¬ 
dence. Hence 

TOc(me) = 1.304(27) GeV. (4) 

is adopted as the final result [ || . In principle the 
same analysis can be performed using the pole 
mass sheme for the quarks. However, the final 
results are quite unstable in contrast to the MS 
scheme. 

The same approach is also applicable to the de¬ 
termination of mb- Again a significant improve¬ 
ment of the stability of the prediction after in¬ 
clusion of the NNLO terms is observed. As final 
result one finds 

mb(mb) = 4.191(51) GeV, (5) 

3. The T(IS') system and the bottom quark 
mass 

In contrast to the previous section one has to 
deal with a non-relativistic system of a bound 
state of a heavy quark-antiquark pair which is 
governed by a complicated multiscale dynamics. 
In the nonrelativistic regime, where the heavy- 
quark velocity v is of the order of the strong¬ 
coupling constant a^, the Coulomb effects are 
crucial and have to be taken into account to all 
orders in Og. This makes the use of the effective 
theory mandatory. This approach allows us to 
separate the scales and to implement the expan¬ 
sion in V at the level of the Lagrangian. 

The dynamics of a nonrelativistic quark- 
antiquark pair is characterized by four different 
regions, the hard region, the soft region, the po¬ 
tential region, and the ultrasoft region. Nonrel¬ 
ativistic QCD (NRQCD) [ ^ is obtained by in¬ 
tegrating out the hard modes. Subsequently in¬ 
tegrating out the soft modes and the potential 
gluons results in the effective theory of potential 
NRQCD (pNRQCD) [ which contains poten¬ 
tial heavy quarks and ultrasoft gluons, ghosts, 
and light quarks as active particles. The effect of 
the modes that have been integrated out is two¬ 
fold: higher-dimensional operators appear in the 


effective Hamiltonian, corresponding to an expan¬ 
sion in V, and the Wilson coefficients of the op¬ 
erators in the effective Hamiltonian acquire cor¬ 
rections, which are series in a^. In [ S I the 
ingredients for the N^LO Hamiltonian have been 
completed using an efficient combination of the 
effective theory formalism and the threshold ex¬ 
pansion [ ^ . 

Once the Hamiltonian is available the only 
task is to solve the corresponding Schrodinger 
equation up to third order using the usual time- 
independent perturbation theory. A new feature 
which appears for the first time at N^LO are the 
retardation effects which arise from the chromo¬ 
electric dipole interaction of the heavy quarko- 
nium with a virtual ultrasoft gluon. 

The third order correction to the lowest energy 
level was calculated in [ P, and will be used in 
this work in order to determine the bottom quark 
mass. The application to the top quark system 
can be found in [ 0 - Recently the N^LO Hamil¬ 
tonian has been used to compute the corrections 
of order In as to the wave function [ . 

The perturbative expansion of the energy level 
with quantum number n looks as follows 

AP-*- =eS + + SEi^'f + 5E^^) + ... , (6) 

where the O (a^) correction, En \ arises from the 
following sources: 

(i) matrix elements of the N^LO operators of 
the effective Hamiltonian between Goulomb 
wave functions; 

(ii) higher iterations of the NLO and NNLO op¬ 
erators of the effective Hamiltonian in time- 
independent perturbation theory; 

(iii) matrix elements of the N^LO instantaneous 
operators generated by the emission and ab¬ 
sorption of ultrasoft gluons; and 

(iv) the retarded ultrasoft contribution. 

(S') 

The analytical result for E)i can be found in [ 
[T^ . In numerical form it reads (adopting the 
choice /is = CFas{fJ,s)'mq) for the bottom and 
top system 

= aUfis)E^ 


f 70.590|,„«„„ 
V 56.7321, 
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Figure 1. mc{mc) for n = 1,2,3 and 4. For each value of n the results from left to right correspond to 
the inclusion of terms of order a°, a\ and in the coefficients Note, that for n = 3 and n = 4 the 
errors can not be determined with the help of Eq. i) in those cases where only the two-loop corrections 
of order as are included into the coefficients Cn as the equation cannot be solved for TOc- 


-1-15.297ln(Q;s(p.s)) -I- 0.001 03 


( 34.229|botto„, 
V 26.654|,„p 



(7) 


where we have separated the contributions arising 
from 03 and / 3 q . The only unknown ingredient in 
the result for 6e[^'^ is the three-loop MS coeffi¬ 
cient 03 of the corrections to the static potential. 
Up to now there are only estimates based on Fade 
approximation [|^ which is used in our analysis. 
However, the final result only changes marginally 
even for a rather large deviations of 03 from its 
Fade estimate. 

The starting point for the determination of the 
bottom quark mass is its relation to the mass of 
the T(IS') resonance 

Mt(is) = 2mf, + Ef *' + , (8) 


with Mr(is) = 9.46030(26) GeV. Here S’^-p Ei 


is the nonperturbative correction to the ground 
state energy. The dominant contribution is asso¬ 
ciated with the gluon condensate and gives [ 0 

P = 60 ± 30 MeV. (9) 

Combining Eq. (H) with the result for the per¬ 
turbative perturbative ground state energy up to 
0{mqal) one obtains the bottom quark mass as 
a function the renormalization scale of the strong 
coupling constant normalization, /r, which is plot¬ 
ted in Fig. 0(a). For the numerical evaluation we 
extract as{mh) with rrih = 4.83 GeV from its 
value at Mz using four-loop /3-fimction accom¬ 
panied with three-loop matchinga. (jUb) is 
used as starting point in order to evaluate (/i) 
at N^LO with the help of the fc-loop f3 function. 

^We use the package RunDec [ to perform the running 
and matching of as 
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(a) 



(b) 


Figure 2. (a) Bottom quark pole mass, mb, as a 
function of the renormalization scale /i, used for 
the numerical evaluation of Ei where the short- 
dashed, long-dashed and solid line corresponds to 
the NLO, NNLO and N^LO approximations, (b) 
MS bottom quark mass rnb(jnb), as a function of 
the renormalization scale n, which is used for the 
extraction of the pole mass mb (cf. (a)). For 
the conversion from the pole to the MS mass the 
method described around Eq. is used. 


From Fig. ||(a) we see that the dependence on the 
renormalization scale becomes very strong below 
/i ~ 2 GeV which indicates that the perturba¬ 
tive corrections are not under control. However, 
even above this scale the perturbative series for 
the pole mass shows no sign of convergence. This 
means that one can assign a numerical value to 


the pole mass only in a specified order of pertur¬ 
bation theory. 

On the contrary, it is widely believed that the 
MS mass mtifJ-) at the scale fi = mb{fi) is a short- 
distance object which has much better perturba¬ 
tive properties. Thus, it seems to be reasonable to 
convert our result for the pole mass into mbirfib). 
The relation between mb and mbirrib) is known up 
to three-loop approximation [ [^, and shows 
sizable perturbative corrections. For this reason 
we suggest the following procedure to take into 
account these corrections in a most accurate way. 
The idea is that from the one-parametric family 
of TOb(/i) we can choose a representative corre¬ 
sponding to some scale /r* in such a way that 

mb{^i*) = mb. ( 10 ) 


For a given fixed-order value of the pole mass 
Eq. (@) can be solved for /r*. In particular, for 
a pole mass to N^LO we use the fc-loop relation 
between the MS and pole mass in Eq. (^. Af¬ 
terwards mbijrib) can be computed from mb{^J.*) 
solving the renormalization group (RG) equation. 
The advantage of this approach is obvious: we 
use the finite order relation between MS and pole 
mass at the scale where they are perturbatively 
close while the large difference between rfibimb) 
and mb is completely covered by the RG evolu¬ 
tion which can be computed with very high accu¬ 
racy as the corresponding anomalous dimension is 
known to four-loop approximation [ ^ M- The 
only restriction on the method could be connected 
to the value of fi*. It should be large enough to 
allow for a reliable use of the RG equation which 
in practice indeed is the case [ 0 - 

In Fig. ib) our result for mb{mb) is plotted 
at NLO, NNLO and N^LO as a function of the 
normalization scale /r which is used to obtain the 
pole mass mb (cf. Fig. §(a)). It is remarkable that 
close to fjL = 2.7 GeV, which is consistent with 
the physically motivated soft scale ~ 2, GeV, 
both the second and the third order corrections 
vanish. This fact is a rather strong indication of 
the convergence of the series for mbirrib). 

The uncertainties in the obtained value of 
mb{rrib) have been investigated in detail in [ 0- 
Finally, the prediction for the MS bottom quark 
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mass reads 

TOfe(m6) = 4.346 ±0.070 GeV. (11) 

4. Conclusions 

For the comparison of the results discussed in 
this contribution with the literature we refer to [ 
I 0- It is, however, interesting to compare 
the value for the charm quark mass with a very 
recent result obtained in a lattice calculation [ 
|T^ . Their final result in quenched approxima¬ 
tion, mc(mc) = 1.301(34) GeV, is impressively 
close to ours (cf. Eq. (||)) with comparable er¬ 
rors. In order to estimate the uncertainty induced 
by the quenched approximation we performed a 
“perturbative quenching” by setting the number 
of active flavours to zero in the calculation of [ Q. 
This leads to a small uncertainty of « 30 MeV, 
which agrees with the estimate of [ 0 - _ 

Although formally slightly beyond la the MS 
quark mass obtained from the low-moment sum 
rule approach is in very good agreement with the 
one from the T(IS') system. One needs to have 
in mind that both the experimental data and the 
theoretical calculations are completely different. 
Whereas no further improvement in the T(15') 
method can be expected there is significant im¬ 
provement possible in the approach discussed in 
Section In particular, after reducing the ex¬ 
perimental error of i?(s) in the charm and bot¬ 
tom threshold region and the one of the leptonic 
widths of the narrow resonances to roughly 2% 
a reduction of the uncertainty in the charm and 
bottom quark mass to 15 MeV and 30 MeV, re¬ 
spectively, can be expected. More details can be 
found in Ref. [ . 
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